Collective modes in the d-density wave state of the cuprates 
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Using a functional integral formulation, we analyze the collective modes in the d-density wave 
state. Since only discrete symmetries are broken, no massless phase mode is present. The only 
relevant fluctuation is the amplitude fluctuation of the order parameter, which is treated at the 
gaussian level above the long range order. These fluctuations give rise to a massive mode, the 
analog of the U(l) Higgs, or amplitude mode at wave-vector (n,n). Neutrons in inelastic neutron 
scattering experiments couple to this mode via the orbital currents and are likely to produce direct 
measurements of this amplitude boson. 
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Much attention has recently been paid to a discrete 
symmetry broken state, the ci-density wave (DDW), as a 
possible candidate for the anomalous aormal state of the 
underdoped cuprate superconductors an Below the pseu- 
dogap temperature scale, T* , the system is conjectured 
to condense into the DDW state, which is responsible for 
much of the observed physics of the pseudogap phase. 
Further down in temperature, below the superconduct- 
ing transition temperature T c , the DDW and the e?-wave 
superconductor (DSC) are expected to coexist and com- 
pete for the same regions on the Fermi surface. In a range 
of the hole-doping concentration, x, suitable for making 
up the underdoped systems, the competition between the 
two orders consistently explains the unique variation of 
the superconducting order parameter with doping]!! 

The DDW state, under various names and, < 
peared in the literature on many occasionsflt 
More elaborate treatments of some of the properties of 
the system, most relevant to experimentSpOn-tiieJiigh-Tc 
superconductors, were discussed recently.EifE3ll3o How- 
ever, these authors, for the most part, treated the system 
in the framework of Hartree-Fock mean field theory, one 
consequence of which was that, all information about the 
collective modes - and experimental implications thereof 
- remained uninvestigated. In this paper, we intend to 
fill this gap by elucidating the structure of the collec- 
tive modes and pointing out how they can actually be 
detected in experiments. 

The DDW state has plaquette-currents, alternating 
clockwise-anticlockwise in the neighboring plaquettes of 
an underlying square lattice. Thus, it gives rise to tiny 
orbital moments. The relevant order-parameter charac- 
terizing the d-density wave is a paricle-hole spin-singlet 
pair, condensing at an wavevector Q = (tt, tt). The bi- 
partite square lattice band-structure is equivalent under 
the transformation Q — » — JO and this forces the order 
parameter to be imaginary,EI a fact which is ultimately 
responsible for the breaking of the time-reversal symme- 
try and the resulting spontaneous bond-currents. The 
phase of the order parameter being fixed by the specific 
symmetries it breaks, the only natural collective excita- 
tion one expects is the amplitude fluctuation of the order- 
parameter, which will give rise to a finite frequency mode, 



the amplitude mode, centered at Q. We will explicitly 
calculate the frequency to be 2Addw, where Addw is 
the DDW single particle gap at (tt,0), which can be es- 
timated from photo-emission experiments above T c . We 
will also show how inelastic neutron scattering, coupling 
to the bond-currents, should directly reveal this mode. In 
fact, well-defined, but damped features at wave-vector Q 
and energy comparable to 2Addw have indeed been seen, 
above T c in inelastic neutron scattering of the cuprates.llj 
These damped peaks are in sharp contrast to the resonant 
peak below T c , which occurs at the same wave- vector and 
comparable energy, but is resolution- limited in frequency. 
This opens up a tantalizing possibility that the amplitude 
mode of the DDW state, which is the analog of the Higgs 
mode of U(l) gauge theory, may actually have already 
been observed in experiments. Clearly, the issue needs 
further investigation. 

In order to analyze the collective modes, let us start 
with the Hamiltonian in momentum space, 
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where creates a particle with momentum k and spin 
a, ek is the free electron band structure —2t(cosk x + 
cos k y ), [X is the chemical potential, and g^ w , is the inter- 
action matrix element. A more realistic band structure 
would include the effect of a next-neighbor hopping if , A 
small t' does not change our results qualitatively, so we 
ignore it for simplicity. In the half filled limit, /i = 0; it 
takes non-zero negative values when holes are doped into 
the system. In this paper we will work in units where the 
lattice constant is unity and will also set H = c = ks = 1. 

The particle-hole spin-singlet d-density wave order pa- 
rameter is anisotropic in fc-space, 
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where $q is the magnitude of the order parameter, and 
/k = {cosk x — cos k y ). So, we will take the effective part 
of the full interaction factorizable as 
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This would enable us to always maintain the d x 2_ y 2-w&ve 
symmetry of the order parameter and discuss long wave- 
length fluctuations around it. For a discussion of the 
DDW state, we will ignore the other parts of the inter- 
action and the fluctuations of other symmetries. 
The Hamiltonian now takes the form 

h = 5> k - aOcLc^ - g- 1 E ^(q)A(q), (4) 



where, 



Note, for q = rQ, a commensurate wave- vector, A(Q) is 
antihermitian,EI and is, in operator form, the DDW gap 
parameter. 

In functional integral form, the partition function is 
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where 



So = J* dr [^{dr + 6 k - /i)4 a c kCT - g- 1 E At(q)A(q)^ , 
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and J\f is a normalization constant. By introducing where ^ is a two-component vector, ^ 
an auxilliary boson field A q through the Hubbard- c£ +Q ), and M is a 2x2 matrix given by, 

Stratonovich transformation, the interaction term is 
made quadratic in the Fermion fields. In frequency space, 
the new action is given by 
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where the cr's are the Pauli matrices acting in the space 
of , 5 kw , and, — (Ji ± i(T2. 

In Eq. ||, the full Brillouin zone has been folded to the 
magnetic, or reduced Brillouin zone (RBZ), to facilitate 
the discussion of the DDW state; the momentum inte- 
grations are restricted to the RBZ, and a spin-sum over 
the Fermionic variables is implied. After integrating out 
the Fermion fields, the partition function is given by, 



z=w f n^dAk,, 
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where, 



S= V-A t (k,w)A(k,w)-21ndetM, (11) 
Q 

and J\f' is now a different normalization constant. The 
factor of two in the second term of S comes from the 



spin-summation. By taking 8S/dA'(\i,u)) — and 
dS/dA(k,uj) = 0, we find the saddle point, 



A(k,w) = 0, fork ^ Q, or lu ^ 0, (12) 



and, 



A(Q,0) = A Q = 5 E/k(^W + £k). (13) 
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Here, -E k = \/ e j + ./u|Aq| 2 and Q(x) is a unit step func- 
tion. Equation |l3j is the d-density wave mean-field gap 
equation. Notice, with this definition of -E k , the DDW 
single-particle gap at (tt, 0), A DDW = 2|Aq|. 

As A(Q,0) is antihcrmitian, the gap parameter 
A(Q,0) must be pure imaginary. Hence, though 
A(Q,0) = Aqe' e satisfies the gap equation equaUjy well, 
the phase of the order parameter is locked to ±§fl Thus, 
there are no massless phase modes in the long- wavelength 
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collective fluctuations of the order parameter. Put an- 
other way, since the order parameter breaks only discrete 
symmetries, there are no Goldstone modes in the system. 
For an incommensurate version of the DDW state, where 
the ordering wave-vector is not commensurate with the 
lattice, there is a collective massless mode, the sliding 
density wave mode. But this mode will be pinned by 
the impurities as well, in second order of the impurity 
potential. 



For an analysis of the amplitude mode, the task at 
hand is to expand the action functional S around the d- 
density wave saddle point. By doing so, we will take 
account of fluctuations of the order parameter ampli- 
tude in the plaquettes, implying fluctuations in the bond- 
currents. This will impose fluctuating non-zero values for 
A(k, oS), at wave- vectors k ^ Q, and frequencies u> 7^ 0. 
Expanding S around the saddle point, to the second or- 
der in A k+Q ^, 
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where the prime over the summation indicates exclusion of the points (k = 0) and (u> — 0). A and B are given by, 
I 



.4 



k+Q,w 



B] 



k+Q,w 



2.9 



k' u>> 



E 



ft' 



(iuo' — — e k /)[i(w' + U)) — fi-\- 



Ck' 



./k'/k'+k 



A 2 
Q 



(15) 
(16) 



Noting that Ak+Q. w is the Fourier component of a pure 
imaginary field, we have the property, 



A_ k+ Q _ u - -Ak + Q.u;. 



(17) 



Using this, and also the fact that A q = — |Aq| 2 in the 
definition of £>k+Q.w, we collapse Eq. |14| onto, 
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where Bk+Q,w is Sk+Q,w with Aq replaced by |Aq| in 
the numerator. 

The dispersion for the collective mode pan be read off 
by equating the fluctuaion kernel to zeroJiS 



to a superconductor, this splitts up into two equations, 
A'(k,u>) - B'(k,u) = 0, and A'(k,u;) + B'(k,u) = 0. 
These are precisely the two equations giving the phase 
and the amplitude mode dispersions, respectively. The 
phase mode is massless, and the massive amplitude mode 
is centered at q = (0,0). In the case of DDW, because 
of the symmetry property of the fluctuation fields which 
leaves the order parameter purely imaginary, one of the 
equations drops out, and, we are left with just one equa- 
tion, that for the amplitude mode. 

Coming back to Eq. |l9L in the case of DDW, writing 
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following the gap equation, in the limit k = 0, the fre- 
quency lo for the amplitude mode is given by, 
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In the case of a superconductor, there is no symme- 
try property, such as Eq. [ij], for the fluctuation fields. 
Gaussian fluctuations around the superconducting saddle 
point has a form very similar to Eq. 14 113 and the col- 



lective modes are given by the determinantal equation, 

= 0. Using the properties, 
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, and B k u = B'_ k _ (i) ,EZl which are specific 
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After the internal frequency summation, the last equa- 
tion, at zero temperature, reduces to, 
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Let us first analyze the mode-frequencies for the sim- 
pler case of an s-wave density wave, the charge density 
wave (CDW). In this case, /k = 1, and Aq = A, a 
real quantity. Even if the order parameter is real, the 
last equation, with the substitution of /k by 1, still gives 
the frequencies of the amplitude mode. This is because, 
for a real order parameter, Eq. holds with a + sign, 
and, -Bk+Q,u and i?k+o,w are same. Consequently, one 
can still collapse Eq. |lj on to Eq. |l8|, and get the same 
equation for the amplitude mode. Considering, then, 
Eq. ^2] for CDW, we note that, for /1 = 0, the roots 
are at lu = ±i2A. This implies the mode- frequencies 
at E = iuj = ±2A. The modes come in pair, for the 
absorptive and emissive responses, respectively. For in- 
creasing values of \fj,\, the modes are shifted, finally, for 
large values of the frequencies scale with ±2 v / 2|/i|- 

Returning to the case of DDW, where /k is nontrivial, 
we have to perform the internal momentum integral nu- 
merically. First, after some simple manipulations, Eq. |2^ 
is brought to the dimcnsinlcss form, 
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where, z = 
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In analogy with the s-wave case, one can easily see 
from here that, at large /u/|Aq|, where the term in 
the numerator can be neglected in comparison, the roots 
of the last equation are at z ~ ±i\/2 ^^ , i.e., at real fre- 
quencies E = iu ~ ±2v2~£t. If n behaves with doping of 
holes in the system in a manner similar to a Fermi liquid, 
then it can be quite large, and eventually scales with t, 
the largest scale in the problem. In such a case, the am- 
plitude mode in the DDW system is at very high energy 
and is virtually inaccesible in thermal neutron scatter- 
ing. Also, a mode so high up in energy is expected to 
be heavily damped, since there is a large phase-space 
of excitations it can decay into. Consequently, one can 
not in the case of large \u\ meaningfully talk about the 
amplitude mode in connection with experiments. But 
a large value of the chemical potential may not-.be the 
case in the underdoped regime of the cuprates JI3 where 
various kinds of charge inhomogeneities and impurities 
may actually pin the doped holes. Indeed, there are some 
indicationsEa that /i remains pinned close to zero in these 
systems. Under such a scenario, the shift of the mode- 
frequencies due to chemical potential is small, and the 
results are qualitatively similar to the case of half-filling, 
u = 0. Thus, in order to make contact with experiments, 
we will only discuss the case of chemical potential pinned 
to zero, assuming it does not shift appreciably with dop- 
ing in the underdoped regime of the cuprates. 

Assuming ji = 0, and replacing z by iu in the left 
hand side of Eq. |2^, we have numerically performed the 
momentum integral over the RBZ. The reason for replac- 
ing z by iu is that, we will eventually look for roots of 
the equation on the imaginary axis such that the mode- 



frequencies after analytic continuation become real. In 
figure 1, we show the results of the integration plotted as 
a function of u. 
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FIG. 1: Results of the integral given in Eq. [23|, I(u), plotted 
as a function of u, where u = —iz. The parameters used are 
t = 0.25 eV, |Aq| = 0.01 eV, and u = 0. I(u) has a root at 
u — 2. Since it is an even function of u, the other root, not 
shown here, occurs at u — — 2. 

Calling the dimensionless function of u in Eq. p3| I{u), 
we note that it has a root at u = 2. Since I(u) is an even 
function of u, there is another root at u = —2. Tracing 
back the transformations, we find the amplitude mode 
frequencies at 
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i2\A Q \z = -2|A Q |« = ±4|A Q | = ±2|A DDW | 



The last equality follows by recalling the DDW single- 
particle gap at (7r,0), Addw = 2| Aq|. So, we conclude 
that the amplitude modes for the DDW state lie at fre- 
quencies ±2Addw- 

At the level of our computation, the mode is un- 
damped, i.e., an exact eigenstate of the system. But, 
since the ordered state is a particle-hole condensate in 
contrast to a superconductor, it can easily couple to other 
degrees of freedom in the cuprates. The mode can, for 
example, couple to the optical phonons, which are abun- 
dant in the system. These other degrees of freedom will 
act as an external bath to the system of the amplitude 
boson, and provide it with some frequency width. We 
can meaningfully talk about the mode - and its manifes- 
tations in experiments - only as long as this width is not 
too big. 

In order to make contact with the neutron scattering 
experiments, let us start with the interaction potential 
between the incoming neutron and the electrons carrying 
the orbital current, V^„tJl2l 
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Vint = ^ tc 

m 
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dl) + h.c. 



(24) 



Here, A is the electromagnetic gauge potential gener- 
ated by the neutron magnetic moment /i n , A = ji n x (r e — 
r n )/(\r e -r n \ 3 ), fj, n = -1.91^S, and, r e and r„ are the 
electron and neutron coordinates, respectively. Because 
the gauge field generated by the neutron is weak, we can 

approximate V int by, V lnt = ietJ2{ij){Ji A ' dl )( c L c j^ ~ 

To find the cross section, we first have to evaluate the 
matrix element of the interaction potential between the 
initial and final states of the neutron, (k^|T^ n t|k,*) = 
/ d Z T n V int expzq • r„, where, q = k 4 - k/._After do- 
ing the integral over the neutron coordinates, til summing 
over the initial and final states in the scattering process, 
and then, putting q = Q = (tt, n), we find, 
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(25) 

Apart from some coefficients, the expression on the 
right hand side of the last equation is precisely the DDW 
gap in the operator form. In other words, the coupling 
of the neutron with the system at (ir, ir) is, as expected, 
through the orbital currents. 

Finally, averaging over the neutron spin states using 
J2 ai a f P^i( a i\^Wf)(gj\^f3\(Ti) = ^ 2 <5 Q/3 , and then, foil- 
wing standard steps,E3 we derive, 



analysis of the DDW amplitude mode. The quantity 
(Aq !W A_q._ w ), in imaginary time ordered form, can be 
directly taken over from Eq. [jj, 
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As we found earlier, this function has poles at u) = 
±i2Addw- These are the dissipative and the absorp- 
tive responses of the system, respectively. At zero tem- 
perature, the relevant pole is that for the dissipative re- 
sponse only. After analytically continuing the frequency, 
the corresponding retarded function will have a pole at 
lu = 2Addw, implying a sharp peak at that frequency. 
This peak captures the amplitude boson for the DDW 
state in inelastic neutron scattering. As we emphasized 
before, the sharpness of the response is really an arti- 
fact of our calculation, which excludes the coupling of 
the boson to other degrees of freedom, such as the op- 
tic phonons. A fuller treatment should bring in some 
damping to the mode. Indeed, the peaks abov£_jT c in the 
neutron scattering of the cuprates are broad.tij Clearly, 
more work is needed in this direction to establish this. To 
establish the correct momentum dispersion will require 
a more microscopic model than that considered. From 
translational invariance one can argue that the square of 
the mode frequency is an analytic function of the mo- 
mentum, that is, it must begin quadratically. However, 
the coefficient of the leading momentum dependence is 
expectedlij to be complex implying damping of the col- 
lective mode at a finite wavevector away from Q = (tt, 7r). 
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Here, Aq iUJ is the frequency transform of the gap pa- 
rameter, Aq, given, in operator form for a general q, in 
Eq. § 

The expression on the right hand side of the last 
equation is precisely the quantity we evaluated in the 
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